AD-A243  501 


ELrCT!^ 

DECi  7  1991 


c 


NOTE 

The  Relationship  Between  Numerical  and 
Physical  Models  of  Atmospheric  Flow 


(Submitted  to  Environmental  Software) 


Roger  A.  Pielke 
Colorado  State  University 
Fort  Collins,  CO  80523 


p 


1  Introduction^ 

There  are  two  fundamental  methods  of  simulating  atmospheric  flows 
-  physical  models  and  mathematical  models.  With  the  first  technique, 
scale  model  repUcas  of  observed  ground  surface  characteristics  (e.g.,  to¬ 
pographic  relief,  buildings)  are  constructed  and  inserted  into  a  chamber 
such  as  a  wind  tunnel  or  a  water  tank.  The  flow  of  air  or  other  gases 
or  liquids  in  this  chamber  is  adjusted  so  as  to  best  represent  the  larger 
scale,  observed  atmospheric  conditions.  Mathematical  modeling,  by  con¬ 
trast,  utilizes  such  basic  analysis  techniques  as  algebra  and  calculus  to 
solve  directly  the  conservation  laws  of  motion,  heat,  moisture,  and  other 
atmospheric  constituents. 

2  Physical  models 

Using  order  of  magnitude  estimates  for  the  dependent  variables  and  as¬ 
suming  that  L  and  S  are  the  representative  length  and  velocity  scale  of 
the  circulation  of  interest  (i.e.,  ^  =  maximum  of 
then  a  scaled  version  of  the  conservation  of  motion  equation  can  be  writ¬ 
ten  as 


where  a  circumflex  (  * )  over  a  dependent  or  independent  variable  indi¬ 
cates  that  it  is  nondimensional.  The  scaling  parameter  e^,  is  a  measure 
of  the  subgrid  scale  velocity  correlations  that  can  be  estimated  from  the 
mean  subgrid  scale  kinetic  energy  i.e.. 

Including  an  estimate  for  the  molecular  viscous  dissipation  and  mul¬ 
tiplying  (1)  by  L/S^  (to  obtain  a  nondimensional  equation  for  the  local 
acceleration)  results  in 

+  (2) 
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To  use  a  scaled  physical  model  to  represent  accurately  the  conservation- 
of-motion  relation  in  the  atmosphere,  it  is  essential  that 

1.  the  individual  bracketed  terms  be  equal  in  the  model  and  in  the 
atmosphere,  or 

2.  the  bracketed  terms  that  are  not  equal  must  be  much  less  in  mag¬ 
nitude  than  the  other  bracketed  terms  in  Elq.  (2). 


I  When  these  conditions  are  met  the  actual  and  modeled  atmospheres  are 
said  to  have  dynamical  similarity. 

Two  of  these  bracketed  terms  are  defined  as 

CIL/S  =  \/Rc 

and 

v/LS  =  1/Re, 

where  Rc  is  the  Rossby  number  and  Re  the  Reynolds  number;  and 

gL(Se/S„)/S^  =  mb 

is  the  bulk  Richardson  number.  (66  represents  the  potential  temperature 
perturbation  and  is  the  same  order  as  the  temperature  perturbation  6T). 

From  Eq.  (2),  to  maintain  dynamic  similarity,  it  is  impUed  that  to 
represent  all  of  the  terms  in  the  equation  properly: 

1.  the  ratio  of  the  subgrid  scale  kinetic  energy  to  the  grid-volume 
average  kinetic  energy  must  be  kept  constant; 

2.  reducing  the  length  scale  L  in  the  physical  model  reqnires; 

(a)  an  increase  in  the  magnitude  of  the  horizontal  temperature 
perturbation  69  or  a  reduction  in  the  simulated  wind  flow 
speed  S  or  both, 

(b)  an  increase  in  the  rotation  rate  fl  or  a  reduction  in  S  or  both, 

(c)  a  decrease  in  the  viscosity  i/  or  an  increase  in  S  or  both; 

3.  an  increase  in  69  in  the  pressure  gradient  term  necessitates  that  S 
also  increase. 

Unfortunately,  it  is  impossible  to  satisfy  all  of  these  requirement  si¬ 
multaneously  in  existing  physical  models  of  mesoscale  atmospheric  circu¬ 
lations.  Such  physical  models  are  constructed  inside  of  buildings,  which 
limits  the  dimensions  of  the  simulated  circulations  to  the  size  of  meters, 
whereas  actual  mesoscale  circulations  extend  over  kilometers. 

To  illustrate  the  difiiculty  of  obtaining  dynamic  similarity  in  a  phys¬ 
ical  model  for  all  terms  in  Eq.  (2),  let  the  horizontal  scale  of  a  mountain 
ridge  be  10  km,  whereas  the  physical  model  of  this  geographic  feature  uti¬ 
lizes  aim  representation.  The  scale  reduction  is,  therefore,  10*.  Thus, 
if  5  =  10  m  s~*  in  the  real  situation  and  air  is  used  in  the  scaled  model 
atmosphere,  then  the  simulated  wind  speed  would  have  to  be  10^  m  s~’ 
to  maintain  identical  Reynolds  number  similarity!  In  addition,  to  have 
the  sanae  Rossby  number  for  this  example,  the  physical  model  must  ro¬ 
tate  10,000  times  more  rapidly  than  the  earth  or  the  wind  speed  must 
be  reduced  by  10,000.  Reducing  the  speed,  of  course,  is  contradictory 
to  what  is  required  to  obtain  Reynolds  number  similarity'  Only  if  the 
results  are  relatively  insensitive  to  changes  in  these  nondimensional  quan¬ 
tities,  as  suggested,  for  example,  by  Cermak  (1975)  for  large  values  of 
the  Reynolds  number  in  simulations  of  the  atmospheric  boundary  layer, 
can  one  ignore  large  differences  in  the  nondimensional  parameters. 

1 

’This  note  wu  extracted  asd  modified  somewhat  trom  that  preaeated  ia  Pielke 
(19M,  Chapter  5). 
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From  the  example  just  given,  however,  it  should  be  clear  that  it  is 
impossible  to  obtain  exact  dynamic  similarity  between  mesoscale  atmo¬ 
spheric  features  and  the  physical  model  when  all  of  the  terms  in  Eiq.  (2) 
are  included. 

Using  the  same  assumptions  applied  to  produce  Eiq.  (1),  the  conserva- 
tion-of-heat  relation,  represented  by  the  potential  temperature  equation, 
can  be  written  as 


where  e»e„,  is  a  measure  of  the  subgrid  scale  correlation  between  the 
fluctuating  velocities  and  temperatures,  with  e*  perhaps  represented  by 


e*  = 

If  the  molecular  conduction  of  potential  temperature  Ct  is  included  in 
(3)  and  represented  in  analogy  with  the  viscous  dissipation  term  as 


C,= 


k,  d^0 

pCj,  dxjdxj ' 


IC.I 


kt  66 


where  kg  is  the  potential  temperature  molecular  conduction  coefficient, 
then  multiplying  by  L/SSS  and  including  the  order  of  magnitude  esti¬ 
mate  of  C$  yields 


temperature  profiles  scaled  according  to  the  nondimensional  relations 
given  by  Eqs.  (2)  and  (4)  and  that  the  flow  will  be  close  to  equilibrium 
(i.e.,  dtildi  and  d6/di  are  small,  relative  to  the  remaining  terms  in 
Eqs.  (2)  and  (4)).  In  addition,  such  bottom  conditions  as  surface  tem¬ 
perature  and  aerodynamic  roughness  must  be  scaled  so  as  to  produce 
kinematic,  dynamic,  and  thermal  similarity  in  the  lowest  levels  of  the 
physical  model.  These  requirements  are  referred  to  as  boundary  ximilar- 
ity  and  their  creation  necessitates  a  comparatively  long  fetch  from  the 
input  region  of  the  laboratory  apparatus  to  the  region  of  simulation,  as 
well  as  obstacles  such  as  a  lattice  placed  upwind  in  the  flow  to  generate 
specific  velocity  profiles  and  turbulence  characteristics. 

With  all  of  these  requirements,  physical  modeling  of  the  mesoscale 
has  been  primarily  limited  to  stably  stratified  flows  over  irregtilar  ter¬ 
rain.  Even  for  this  case,  however,  such  observed  features  of  the  real 
atmosphere  as  the  veering  of  the  winds  with  height,  radiational  cooling, 
and  condensation  cannot  be  reproduced. 

3  Use  of  numerical  models  to  provide  bound¬ 
ary  conditions  for  the  physical  models 

Pielke  (1984)  reviews  in  depth  the  techniques  of  numerical  modeling  of 
atmospheric  flows.  For  the  purpose  of  this  paper,  such  numerical  models 
can  be  used  to  remove  some  of  the  constraints  placed  on  physical  models 
by  providing; 


,  39 

di  ~ 


dxj  ^ 


(4) 


The  ratio 


ki/pdCgO  =  Pr-‘ 

is  called  the  Prandil  number  and  is  of  order  unity  for  air. 

Thus  to  obtain  thermal  similarity  between  the  mesoscale  drculation 
and  its  laboratory  representation,  the  Reynolds  number  must  also  be 
very  large  and  the  partitioning  of  heat  transport  between  the  subgrid 
scale  and  resolvable  fluxes  must  be  the  same.  If  in  Eq.  (2),  for  ex¬ 
ample,  the  temperature  perturbation  66  must  be  increased  in  the  bulk 
Richardson  number  Ribuik  to  compensate  for  a  decrease  of  L  in  the  labo¬ 
ratory  model,  then  in  Eq.  (4)  the  turbulent  fluctuations  in  the  simulated 
atmosphere  must  also  be  increased. 

The  nondimensional  source-sink  term  for  potential  temperature  Ss 
is  included  in  the  analysis.  However,  the  mathematical  procedure  of 
representing  it  as  a  single  variable  masks  its  physical  complexity.  This 
term  includes  such  effects  as  radiative  flux  divergence,  phase  changes  of 
water,  etc.  and  is  an  involved  function  of  the  dependent  variables.  Thus 
it  is  extremely  difficult  to  evaluate  this  term  using  scale  analysis,  and, 
in  practice,  physical  modelers  exclude  it  in  their  representation  of  at¬ 
mospheric  flows.  An  equivalent  similarity  analysis  can  be  performed  for 
water  substance  and  other  aerosol  and  gaseous  contaminants.  Because  of 
the  inability^  to  accurately  represent  the  sources  and  sinks  of  these  vari¬ 
ables,  (i.e.,  Sq,  and  5i„),  however,  physical  modelers  have  only  studied 
the  movement  of  nonreactive,  conservative  pollutants  around  terrain  and 
building  obstacles. 

In  utilizing  physical  models,  the  conservation-of-mass  relation  must 
also  be  satisfied,  and,  as  long  as  the  ratio  of  the  variations  of  spe¬ 
cific  volume  to  the  average  specific  volume  in  the  physical  model  are 
much  less  than  unity  explicit  temporal  changes  in  density  can  be  ignored 
in  this  conservation  relation.  The  scaled  version  of  the  incompressible 
conservation-of-mass  equation,  =  0,  shows  that 

yv  ~  L.V/L^,  W  ~  L.V/L, 

io  that,  if  the  ratio  of  the  vertical  to  the  horizontal  scales  of  the  circu- 
ation  is  kept  constant  between  the  physical  model  and  the  atmosphere, 
kinematic  similarity  is  obtained.  This  requirement  could  be  satisfied  pro¬ 
viding  the  horizontal  to  vertical  representation  of  the  terrain  and  other 
>hysical  features  of  the  ground  surface  in  the  physical  model  are  not 
ixaggerated.  This  latter  condition  is  called  geometric  similarity. 

The  final  similarity  conditions  needed  in  physical  models  include  the 
equirement  that  air  flowing  into  the  simulated  region  have  velocity  and 


1) .  the  appropriately  scaled  velocity,  temperature,  moisture,  and  other 

atmospheric  constituent  profiles  in  the  inflow  stream  to  the  physical 
model; 

and 

2) .  the  appropriate  surface  fluxes  of  heat,  moisture  and  momentum,  as 

scaled  for  the  physical  model  (e.g.,  eJ,/5*  from  (2),  and  ese„l66S 
from  (4)). 

In  this  approach,  the  numerical-physical  models  are  linked  in  a  one-way 
nested  grid  approach  where  the  numerical  atmospheric  model  provides 
a  realistic  flow  simulation  for  a  larger  scale  domain  at  a  scale  where  the 
physical  model  suffers  from  similarity  limitations.  This  time-evolving 
flow  informatiem  is  used  as  boundary  conditions  for  the  smaUer  domain 
physical  model.  The  physical  model  can  provide  realistic  simulations  on 
the  smaller  scale  where  the  numerical  model  suffers  from  computational 
limitations  due  to  the  small  domain  size  and  resultant  large  sensitivity 
of  the  results  to  even  slight  errors  in  the  lateral  boundary  conditions  of 
the  numerical  model  (Pielke,  1987).  In  addition,  the  numerical  models 
fail  to  adequately  resolve  nearly  discontinuous  obstacles  such  as  cliffs, 
buildings,  individual  trees,  etc. 
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List  of  Symbols 


Dependent  Variables 

Uj  =  velocity  vector 
8  =  potential  temperature 

T  =  scaled  pressure;  r  =  Cp(p/1000  where  Cp  is  the 

specific  heat  at  constant  pressure,  p  is  pressure  in  millibars 

Constants 

R  =  gas  constant 

n  =  rotation  rate  of  the  earth  in  radians  per  second 
1/  =  kinematic  viscosity  of  air 
g  =  gravitational  acceleration 

kf  =  molecular  coefficient  of  conduction  of  potential  temper¬ 
ature 

Mathematical  Operations 

(")  =  grid-volume  average  of  dependent  variable,  e.g., 
etc. 

(  )'  =  resolvable  perturbation  of  dependent  variable  from 
domain-average  variable 

( )o  =  domain  average 

(*)  =  nondimensional  variable 

(  )"  =  subgrid  scale  perturbation  of  dependent  variable  from 
grid-volume  average 

Scaling  Values 

L  =  length  scale 

68  =  scaled  temperature  gradient 
5  =  scaled  velocity 

Independent  Variable 

Xi  =  spatial  independent  variable 
(  =  time 
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